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An algorithm based on quantum phase estimation, which discriminates quantum states non- 
destructive ly within a set of arbitrary orthogonal states, is described and experimentally verified by 
a NMR quantum information processor. The procedure is scalable and can be applied to any set of 
orthogonal states. Scalability is demonstrated through Mat lab simulation. 



I. INTRODUCTION 

Quantum Computing has generated re-newed interest 
in theory and practice of Quantum Mechanics. This is 
largely due to the fact that quantum computers can solve 
certain problems much faster than classical computers 
[T]-[3j. Many efforts are being made to realize a scal- 
able quantum computer using techniques such as trapped 
ions, optical lattices, diamond-based quantum comput- 
ers, Bose-Einstein condensate based quantum comput- 
ers, cavity quantum electrodynamics (CQED) and nu- 
clear magnetic resonance ^-0. NMR has become an 
important experimental tool for demonstrating quantum 
algorithms, simulating quantum systems, and for verify- 
ing various tenets of quantum mechanics [91420] . 

There are several theoretical protocols available for 
orthogonal state discrimination [2TH24] . Walgate et al. 
showed that, using local operations and classical com- 
munication (LOCC) multipartite orthogonal states can 
be distinguished perfectly [21]. However if only a single 
copy is provided and only LOCC is allowed, it cannot 
discriminate quantum states either deterministically or 
probabilistically [22]. Estimation of the phase plays an 
important role in quantum information processing and 
is a key subroutine of many quantum algorithms. When 
the phase estimation is combined with other quantum al- 
gorithms, it can be employed to perform certain compu- 
tational tasks such as quantum counting, order finding 
and factorization [2^, ^25^. Phase Estimation Algorithm 
has also been utilised in a recent important application 
in which the ground state of the Hydrogen molecule has 
been obtained upto 45 bit accuracy in NMR and upto 20 
bit accuracy in photonic systems [26 . 

By defining an operator with preferred eigen- values, 
phase estimation can be used logically for discrimina- 
tion of quantum states with certainty [6 . It preserves 
the state since local operations on ancilla qubit mea- 
surements do not affect the state. In this paper we de- 
scribe an algorithm for non-destructive state discrimina- 
tion using only phase estimation alone. The algorithm 



described in this paper is scalable and can be used for 
discriminating any set of orthogonal states (entangled or 
non-entangled). Earlier non-destructive Bell state dis- 
crimination has been described by Gupta et al. [27 and 
verified experimentally in our laboratory by Jharana et 
al. [23. Bell states are specific example of orthogonal 
entangled states. The circuit used for Bell state dis- 
crimination [27] is based on parity and phase estimation 
and will not be able to discriminate a superposition state 
which has no definite parity. For example consider a state 
= -^(|00) + 1 01)), which belongs to a set of orthog- 
onal states. Here |00) has parity and |01) has parity 
1. Hence the above lip) does not have a definite par- 
ity and cannot be distinguished from its other members 
of the set, by the method of Gupta et al. [2T. Secjllj 
of this paper describes the design of a circuit for non- 
destructive state discrimination using phase estimation. 
Sec. ini also contains non-destructive discrimination of 
special cases such as Bell states and three qubit GHZ 



states using phase estimation. Sec|III| describes exper- 
imental implementation of the algorithm for two qubit 
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states by NMR quantum computer and Sec IV describes 
the Matlab® simulation of non-destructive discrimina- 
tion of three qubit GHZ states. 



II. THEORY 

For a given eigen- vector of a unitary operator [/, 
phase estimation circuit with Controlled-U operator can 
be used for finding the eigen- value of |0) [6]. Conversely 
the reverse of the algorithm, with defined eigen-values 
can be used for discriminating eigen-vectors. By logi- 
cally defining the operators with preferred eigen-values, 
the discrimination, as shown here, can be done with cer- 
tainty. 



A. The General Procedure (n-qubit case): 

For n qubit case the Hilbert space dimension is 2^, 
having 2^ independent orthogonal states. Hence we need 
to design a quantum circuit for state discrimination for 
a set of 2^ orthogonal quantum states. Consider a set 
of 2^ orthogonal states {^i}, where i = 1,2, ....2^. The 
main aim of the discrimination circuit is to make direct 
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where V is the matrix formed by the column vectors 

v= 1,^2) \^^) !</-„)]. 

The circuit diagram for implementation of Phase Es- 
timation Algorithm (PEA) to discriminate orthogonal 
states using the Controlled-Uj operations such that the 
original state is preserved for further use in any qauntum 
circuit is shown in Fig.l. 

As the eigen- values defined are either +1 or —1, the 
final ancilla qubit states will be in product state (with- 
out superposition) , and hence can be measured with cer- 
tainty. It can be shown that the selection of specific 
operator set {Uj} with the conditions discussed above 
makes direct correlation between 2^ product states of 
ancilla qubit and elements of so that ancilla mea- 

surements can discriminate the state. 



B. Single qubit case: 

For a single qubit system, the Hilbert space dimension 
is 2. So we can discriminate a state from a set of two 
orthogonal states. Consider an illustrative example with 
the orthonormal set as { = ^^^dO) + ll)), \4^2) = 
^(|0) — |1)) }. The quantum circuit for this particular 
case can be designed by following the general procedure 
discussed in Sec II A| The V matrix for the given states 
{\^i),\4>2)} is, 



correlation between the elements of {^i} and possible 
product states of ancilla qubits. As there are 2^ states, 
we need n ancilla qubits for proper discrimination. 

The discrimination circuit requires n Controlled 
Operations. Selecting these n operators {Uj} (where 
j = 1, 2, ...n) is the main task in designing the algorithm. 
The set {Uj} depends on the 2^ orthogonal states in 
such a way that the set of orthogonal vectors forms the 
eigen- vector set of the operators, with eigen- values ±1. 
The sequence of +1 and —1 in the eigen- values should 
be defined in a special way, as outlined below. Let {e*} 
(with i = 1, 2... 2^) be the eigen- value array of Uj^ and it 
should satisfy following conditions. 

Condition #1: Eigen- value arrays {e*} of all operators 
{Uj} should contain equal number of +1 and -1, 
Condition #2: For the first operator /7i, the eigen- value 
array {e^} can be any possible sequence of +1 and -1 
with Condition #1^ 

Condition #3: The restriction on eigen- value arrays 
starts from Uj=2 onwards. The eigen- value array ({62}) 
of operator U2 should not be equal to {e^} or its 
complement, while still satisfying the Condition #1. 
Condition #4 ' By generalizing the Condition #3, the 
eigen- value array ({e^}) of operator Uk should not be 
equal to {e^} (m = 1, 2, ...k — 1) or its complement. 

Let Mj be the diagonal matrix formed by eigen-value 
array {e*} of Uj. The operator Uj is directly related to with 
Mj by a unitary transformation given by. 
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According to the rules given in Sec jll A[ M can be either 



U = V-' X M xV = 



1 

1 



(2) 



The circuit diagram for this case is identical to Fig.l, 
having only one work and one ancilla qubit. It can be 
easily shown that, the ancilla qubit measurements are 
directly correlated with the input states. For the selected 
Ml, if the given state is then ancilla will be in the 
state |0) and if the given state is \(j)2) ancilla will be in 
the state 

For a general set {|^i) = (a|0) + /3|1)), 
1^2) = (/^|0) ~ <^|1))} (where a and /3 are real numbers 
satisfying, + = 1), operator U for eigenvalue 
array {1,-1} can be shown as. 



U-- 



fCos{0) Sin{0) \ 
\Sin{e) -Cos{0)J 



(3) 
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C. Two qubit case: 

The Hilbert space dimension of two qubit system of 
is four. Consider an illustrative example with a set of 
orthogonal states 



{\S{a,m 



{{am 
(/?|10) - 



f/?|01)),(a|10) 
a|ll)),(/?|00)- 



a\01))}, (4) 



where a and /3 are real numbers satisfying, + = 
1. This set is so chosen that the states are (a)orthogonal, 
(b)not entangled, (c)different from Bell states, (d)do not 
have definite parity and (e) contain single-superposed- 
qubits (SSQB) (in this case second qubit is superposed). 
Using the general procedure discussed above, we can se- 
lect the eigen-value arrays for two operators Ui and U2 
as 



{ei} = {l, 1,-1,-1}, 



{62} = {1,-1, 1,-1}. (5) 



Ui and U2, the unitary transformation of the diagonal 
matrices formed by {ei} and {62} are. 



Ui = 



/Cos{0) Sin{0) 

Sin{e) -Cos{e) 

Cos{e) Sin{0) 

V Sin{e) -Cos{0)/ 



(6) 



Ancilla Qubits 




FIG. 1. The general circuit for non-destructive Quantum State Discrimination. For discriminating n qubit states it uses n 
number of ancilla qubits with n controlled operations, n ancilla qubits are first prepared in the state |00...0). Here H represents 
Hadamard transform and the meter represents a measurement of the qubit state. The original state encoded in n qubits is 
preserved (not destroyed). 
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TABLE I. State of ancilla qubits for different input states for 
two qubit orthogonal states. 



U2 = 



fCos{e) 

Sin{0) 


V 



Sin{9) 
-Cos{0) 
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-Cos{e) -Sin{e) 

-Sin{9) Cos{e) j 



(7) 



where, 9 = 2 x Tan ^ ( 



/3^ 



The output state of the ancilla qubit run through all 
possible product states as input state changes, as listed 
in TableJ. The quantum circuit for two qubit state dis- 
crimination is shown in Fig. 2a. 



1. Special case (a = /3 = -^): 

The set of orthogonal states are, 

{|5(^,^))} = {|0,)} = {^(|OO) + |O1)), 
^(|10) + 111)), ^(|10) - 111)), ^(100) - 101))}. (8) 

The operators Ui and U2 can be found by substituting 
the value of ^ = ^ in (5) and (6), 
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and U2 
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(9) 

The quantum circuit for the set (Eqn.8) is same as 
the general case of any set of two qubit orthogonal 
states(Fig.2a). Experimental implementation of this case 
has been performed using NMR and is described in Sec. 

m 



D. Bell state discrimination: 

Bell states are maximally entangled two qubit states 
(also known as Einstein- Podolsky- Rosen states) [29] . 
They play a crucial role in several applications of 
quantum computation and quantum information theory. 
They have been used for teleportation, dense coding and 
entanglement swapping [W-'ST. Bell states have also 
found application in remote state preparation, where a 
known state is prepared in a distant laboratory [ 34] . 
Hence, it is of general interest to distinguish Bell states 
without disturbing them. The complete set of Bell states 
is, 



{m} = {j=,{\020^) 
102I1)- 
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Bell states form an orthogonal set. Hence one can design 
a circuit for Bell state discrimination using only phase 
estimation. The circuit diagram is same as that shown 
in Fig. 2a with different Ui and 1/2- For eigen- value 
arrays {ei} = {1, -1, 1, -1}, {62} = {-1,1,1,-1}. Ui 
and U2 are obtained as. 
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The controlled operators {C — Ui and C — U2) for phase 
estimation which involves 3-qubit operators can be writ- 
ten as the product of 2-qubit operators as. 



C 
C 



U2 



C - NOTl xC- NOT; 



2 5 



(12) 



C-7tIxC- NOT^ xC- not! X C 



Here qubits 1 and 2 are work qubits in which the Bell 
states are encoded and 3 and 4 are the ancilla qubits. 
Here C—NOTj represents C-NOT operation with control 
on i^^ qubit and target on j^^ qubit. The splitting of 
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three qubit operator into two qubit operators is needed 
for the implement ion of C — Ui and C — U2- 

There aheady exists an algorithm for non-destructive 
discrimination of Bell state by Gupta et al. [27], which 
has also been experimentally implemented in NMR by 
Jharana et al. [28 . The circuit of Gupta et al. [27 is 
based on parity and phase measurement and will fail for 
a superposed state which has no definite parity. However 
for non-destructive discrimination of Bell states using the 
present phase estimation algorithm is similiar to Gupta's 
circuit where the parity estimation is replaced by modi- 
fied phase estimation. 

E. GHZ state discrimination: 

GHZ states are maximally entangled multi qubit 
states [35 . GHZ states have been used in several quan- 
tum algorithms such as quantum secret sharing, con- 
trolled dense coding and quantum key distribution [36]- 
[38]. These algorithms make use of entanglement and 
hence it is important to discriminate GHZ states by pre- 
serving their entanglement. 

All n qubit GHZ states form an orthogonal set (with- 
out definite parity). Hence a circuit can be designed for 
discriminating general n qubit GHZ states using only 
phase estimation. 

Consider the case of three qubit GHZ states, which 
are 

= {;^(|000) + 1111)), ^(1000) - 1111)), 
^(1001) + I110)), ^(1001) -|110)), 
^(1010) + I101)), ^(1010) -|101)), 
^(1011) + 1100)), ^(1011) -1100))} (13) 

Here we need three ancilla qubits and have to imple- 
ment three controlled operators for state discrimination. 
Verification of the NMR experiment to discriminate such 
states has been carried out here using Matlab® and the 
parameters of a four qubit NMR system, as described in 
SecHVl 



III. EXPERIMENTAL IMPLEMENTATION BY 
NMR 

A. Non-Destructive Discrimination of two qubit 
orthogonal states: 

Experimental implementation of the quantum state 
discrimination(QSD) algorithm has been performed here 
for 2 qubit case for the orthogonal set (Eqn.7). 

The discrimination circuit diagram shown in Fig. 2a, 
needs a 4 qubit system. As two ancilla qubits are inde- 
pendent of each other, following [28] one can split the ex- 
periment into 2 measurements with a single ancilla qubit 
(Fig.2b and 2c). 



The NMR implementation of the discrimination al- 
gorithm starts with (i) preparation of the pseudo-pure 
state followed by (ii) creation of input state (iii) Quan- 
tum Phase estimation with operators {Uj}. Finally the 
measurement on ancilla qubits yields the result. 

The experiment has been carried out at 300K in 11. 7T 
field in a Bruker AVSOO spectrometer using a triple res- 
onance QXI probe. The system chosen for the imple- 
mentation of the discrimination algorithm is Carbon- 13 
labeled Dibromo-fluoro methane ^^CHFBr2, where ^H, 
^^F and ^^C act as the three qubits[39l gD]. The ^H, 
^^F and ^^C resonance frequencies at this field are 500, 
470 and 125 MHz, respectively. The scalar couplings be- 
tween the spins are: Jhc = 224.5i^z, Jhf = 49.7i^z 
and Jfc = -310.9i^z(Fig.3). 

The NMR Hamiltonian for a three qubit weakly cou- 
pled spin system is [41 , 

3 3 
i=0 i<j=l 

where Ui are the Larmor frequencies and the Jij are 
the scalar couplings. The starting point of any algorithm 
in an NMR quantum information processor is the equi- 
librium density matrix, which under high temperature 
and high field approximation is in a highly mixed state 
represented bv[40j. 

Peq JhI^ + 7c/f + 7F/f 

= 7//(/f + 0.94/f + 0.25/f ). (15) 

There are several methods for creating pseudo pure states 
(PPS) in NMR from equlibrium state [lOl glHll]. We 
have utilized the spatial averaging technique [10] for cre- 
ating pseudo pure states as described in [40^. The spectra 
for equlibrium and |000) PPS are shown in Fig. 3. 

For Phase Estimation algorithm, due to its high sensi- 
tivity, proton spin has been utilized as the ancilla qubit; 
and the two qubit states, to be discriminated, are en- 
coded in carbon and fluorine spins. As the measure- 
ments are performed only on ancilla qubit, we record 
only proton spectra for non-destructive discrimination of 
the state of carbon and fluorine. The state of the ancilla 
qubit can be identified by the relative phase of the spec- 
tra. We set the phase such that a positive peak indicates 
that the proton was initially in state |0). 

Implementation of Controlled -Ui and U2' 

For the set of orthogonal states given in eqn.(9) the Ui 
and U2 are given in Eqn.(8). Let Hi and H2 be the effec- 
tive Hamiltonians for Controlled-Ui and Controlled-U2 
propagators such that, 
Controlled-Ui = exp{iHi), 
Controlled-U2 = exp{iH2). 

where Hi and 7^2, in terms of product operators [IT" are 
obtained as. 
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Ancilla Qubit -2 

|o> — 



{1^0} 



FIG. 2. (a) Two qubit State discrimination circuit for Experimental implementation in three qubit NMR quantum computer, 
(b) and (c) are splitting of the circuit- (a) into two circuits with single ancilla measurements. 
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FIG. 3. The three qubit NMR sample used for experimental implementation. The nuclear spins ^if, ^^F and ^^C are used as 
the three qubits. (a) Equlibrium Spectra of proton, carbon and fluorine, (b) Spectra corresponds to the created |000) pseudo 
pure state. These spectra are obtain by using 90^^ measuring pulse on each spin. 



H2 = (J/ - |/i - 7r/2/3 + 27r/i/2/3). 

Since the various terms in Hi and H2 commute with each 
other, one can write, 

Controlled — Ui= exp{iHi) 

= exp{i^I) X exp{-i^ll)xexp{-i^I^) x exp{mllll), 
Controlled — U2 = exp{iH2) 

= exp{iC-I - ^ll-T^I^ll + 27r/i/2/3)) 

= exp{i^I) X exp{-i^ll) x exp{-inl^ll) 

X exp(i2TTllllll). (16) 

As the decomposed terms commute with each 
other, these propagators can be easily implemented in 
NMR(Fig.4). Single spin operators such as IrjQ , ly are 
implemented using R.F pulses. The Iz operator is im- 
plemented using composite z rotation pulses in NMR 
{{l)-x{l)y{l)x) [13 Hg. Two spin product terms such 
as Illl are implemented using scalar coupling Hamilto- 
nian evolution sandwiched between two {^)y pulse on 
j spin [40 . The three spin product operator terms are 
implemented using cascades of two spin operator evolu- 



tions( T^en^ et al. [47]). 

The experimental results are shown in Fig. 5. Pro- 
ton spectra shows the state of ancilla qubit, which in- 
turn can be used for discrimination of two qubit state 
in carbon and fluorine spins. Positive peaks in Fig. 5 
means ancilla is in qubit state |0) and negative peak 
indicates ancilla qubit is in state Thus spectra 

in Fig. 5 indicates that (i), (ii), (in), {iv) are respectively 
1^1)7 1</>2), |(/>3), |(/>4) (Table. I). To compute fidelity of the 
experiment, complete density matrix tomography has 
been carried out (Fig. 6). 

The experimental results are in agreement with the 
Table. I with an 'average absolute deviation' [48 of 4.0% 
and 'maximum absolute deviation' [48] of 7.2%, providing 
the desired discrimination. 



IV. THREE QUBIT GHZ STATE 
DISCRIMINATION USING Matlab® 
SIMULATION: 

Non-destructive discrimination of the three qubit max- 
imally entangled (GHZ) states using only Phase Esti- 
mation algorithm as described in Secjiljin NMR has also 
been performed using a Matlab® simulation. This simu- 
lation verifies the principle involved but does not include 
any decoherence or pulse imperfection effects. The three 
qubit GHZ states form a set {Gi} given by eqn.(13) can 
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FIG. 4. The pulse sequence for Controlled-Uj operators for two qubit orthogonal states shown in (14) (Here narrow pulses 
indicate (|) pulses and broad pulses indicate tt pulses with the phase given above the pulse). 



{ii) 



(Hi) 



(iv) 



FIG. 5. Ancilla (proton spin) spectra of final state for two qubit state discrimination algorithm for (i) (ii) |02), (Hi) l^s), 

(iv) 104) states. Ai and A2 are results of two measurements on single ancilla(Fig.4b and 4c respectively) qubit- 1 and 2 (here 
it is two experiments with same ancilla qubit). These spectra are obtained with a 90^ measuring pulse on the ancilla(proton) 
qubit at the end of the pulse sequence. 



be re-expressed as, 

{|Gi)} = {|(/.l),|<^l),|(/.3), |<^8)} (17) 

The discrimination of a 3-qubit GHZ state using phase 
estimation requires 3 work qubits and 3 ancilla. We di- 
vide the 6 qubit quantum circuit into three circuits. Each 
circuit has three work qubits and a single ancilla. There 
are several possibilities for eigen- value sets which will sat- 
isfy the sets of conditions discussed in SecJII| Consider 
one such set, 

{ei} = {l,-l,l, -1,1, -1,1,-1}, 

{62} = {-1,1, 1,-1, 1,-1, -1,1}, (18) 

{63} = {-1,1, 1,-1, -1,1, 1,-1}. 

For this eigen- value set (18), the Controlled — Uj opera- 
tors can be written as 

Controlled - Ui = C - NOT^ xC- NOT^ xC- NOT^, 
Controlled -U2 = C -irlxC - NOT^ xC- NOT^ 

X C-NOT^ xC-nl 
Controlled - U^. = C - ttI x C - NOT^ xC - NOT^ 
xC-NOTl xC-irl. (19) 

Splitting of four qubit operator into two qubit opera- 
tors is needed for its experimental implementation (Fig. 7) 
Here 1,2 and 3 are the work qubits, in which the GHZ 
state is encoded and 'a' is the ancilla qubit. The results 
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|0> 
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TABLE II. State of ancilla qubits for different input states of 
Eqn.(13) and (17). 



of ancilla qubit measurements are tabulated in Table. IT 

NMR simulation has been carried out using the pa- 
rameters of a well known 4-qubit system, crotonic acid 
with all carbons labelled by ^^C(Fig.8) [16 . The den- 
sity matrix tomography of the Matlab® experiment for 
a few selected ( 1^1), 1^4) and 1^7)) GHZ states are shown 
in Fig. 10. This confirms that the method of Phase Esti- 
mation discussed in SecHIlcan be used for discrimination 
of GHZ states without destroying them. 



CONCLUSION 

A general scalable method for non-destructive quan- 
tum state discrimination of a set of orthogonal states 
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Initial State After First Experiment After Second Experiment 

(i) |0>(;^(100> + |01») |0)(-i5(|00) + |01))) loX^ijdOO) + 101))) 




(ii) |0>(;)5(|10) + |11») |0>(-i5(|10> + |ll))) |1)(;)5(|10> + 111))) new file 




(vi) |0)(-l5(|00)-|01))) |l)( 1 (|00)-101») |1)( 1(100) -|01))) 




FIG. 6. Density Matrix Tomography of the initial and final states of QSD circuit. First qubit is the ancilla. It is evident 
that the state of 2^^ and 3^^ qubits are preserved, (here 1 |000), 2 |001), 3 |010), 4 |011), 5 |100), 6 |101), 
7^ |110), 8^ |111).) 



using quantum phase estimation algorithm has been 
descibed, and experimently implemented for a two 
qubit case by NMR. As the direct measurements are 
performed only on the ancilla, the discriminated states 
are preserved. The generalization of the algorithm is 



illustrated by discrimination of GHZ states using a 
Matlab® simulation. 
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FIG. 8. The chemical structure, the chemical shifts and spin-spin coupling of a ^^C labelled Crotonic Acid. The four ^^C spins 
act as four qubits [16j . 
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FIG. 9. Matlah® simulation results for GHZ state discrimination. The simulated spectras are shown for three GHZ states 
and 104)- It is evident from final density matrix that the GHZ states are preserved. (Here 1 |0000), 2 |0001), 
3 ^ |0010), 4 ^ |0011), 5 ^ |0100), 6 ^ |0101), 7 ^ |0110), 8 ^ |0111), 9 ^ |1000), 10 ^ |1001), 11 ^ |1010), 12 ^ |1011), 
13 |1100), 14 |1101), 15 |1110), 16 |1111). First qubit is the ancilla) 



